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Abstract 

Cd ' Statistics of electron tunneling in normal tunnel junctions is studied analytically and numerically 

J^ ■ taking into account circuit (environment) effects. Full counting statistics, as well as full statistics 

of voltage and phase have been found for arbitrary times of observation. The theoretical analysis 
j^ I was based on the classical master equation, whereas the numerical simulations employed standard 

Monte-Carlo methods. 
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I. INTRODUCTION 

Shot noise produced by electron current is known nearly 90 yearsi. The interest to this 
phenomenon was revived recently due to great importance of shot noise for modern meso- 
scopic devices^. Shot noise provides valuable information on the charge of carriers, which are 
responsible for the electric current. Now researchers are not content with the parametriza- 
tion of noise using variance or noise temperature and look for a complete picture of charge 
noise given by all cumulants of the probability distribution for the number of electrons, 
which traverse the junction during a fixed period of time. The full set of moments or cu- 
mulants determines the full counting statistics. The theoretical study of the full counting 
statistics was initiated by Levitov and Lezovik^, and recently these studies have been in- 
tensified and brought about important results concerning the effect of environment on the 
statistics (crossover from the voltage to the current bias)^. These studies focused on the 
long-time (low-frequency) limit of the full counting statistics. 

It is difficult to study full counting statistics experimentally, but essential progress has 
been achieved in this direction too. Experimentalists do not have direct access to counting 
statistics itself, i.e. they cannot determine those moments of time when electrons cross the 
junction, but instead they can scan voltage noise produced over a shunt by tunneling events. 
Until now, only the first non-trivial cumulant, namely the third cumulant (skewness)§!fl, has 
been detected experimentally. In the literature, there has been discussion of other methods 
of noise spectroscopy. It has been demonstrated theoretically and experimentally that a 
Coulomb-blockade tunnel junction is an effective probe of shot noise^ (and of other types of 
noise as well^) . On the other hand, Tobiska and Nazarov^ suggested to use a superconducting 
Josephson junction close to the critical current as a threshold detector. This method has 
also been studied experimentally^. 

In tunnel junctions, shot noise is produced by discrete electrons, which tunnel quantum- 
mechanically through a high potential barrier and, thus, noise characterizes this quantum- 
mechanical process. This was the reason to investigate shot noise with modern quantum-field 
techniques^. However, it is known that though quantum mechanics is crucial for formulation 
of basic statistic properties of electron transport, after the basic statistical properties of 
electron tunneling have already been formulated, the following statistical analysis, which 
should lead to knowledge of full statistics can be done without any reference to quantum 



mechanics. In particular, Nagaevii have studied the effect of environment (circuit) on the 
counting statistics using the classic Boltzman-Langevin approach, which agrees with the 
results of the quantum-mechanical analysis. 

Our present work describes a further development of the classical analysis of statistics of 
electron tunneling using the formalism of the master equation and direct numerical simula- 
tions. The master equation has been widely used for studying various problems of statistics 
in physics and in other fields. The master equation for the voltage probability distribution 
was used long time ago within the framework of the "orthodox theory" of the Coulomb 
blockadei^. But it was mostly for the calculation of the dynamics of single electron tun- 
neling oscillations and the corresponding IV curves at low bias, when Coulomb blockade 
modifies the IV curve essentially. The master equation has also been employed for studying 
counting statistics in various mesoscopic setups^^, but only in the long-time (low-frequency) 
limit. Here our intention is to use the master equation for full statistics of charge transfer 
through a normal junction for any time scale. We write here full statistics instead of full 
counting statistics since we have analyzed not only the statistics of electron tunneling events 
(counting) but also the statistics of voltage and phase generated by these events. This is 
important, since, as mentioned above, typical experiments probe voltage or phase instead 
of the number of electrons. In the long-time (low frequency) limit, our results completely 
agree with the previous quantum-mechanical analysis^. However, our final expressions valid 
at T = are not restricted to long times and, therefore, they provide the full statistics for 
all time scales. 

II. ELECTRIC CIRCUIT AND IV CURVE 

We consider the simplest circuit, which is standard for studies of shot noise in tunnel 
junctions: a tunnel junction of resistance Rt and of capacitance C biased from an ideal 
voltage source Vq via a resistance R (called shunt) in series (Fig. [T^). The voltage across 
the junction is Vt = Vq — V where V = IR is the voltage over the shunt resistor and 
/ = Vq/{Rt + R) is the current through the junction. A basic parameter for our analysis 
is the ratio a = R/ Rt- If a <^ 1 the circuit corresponds to an ideal voltage bias Vq ~ Vr, 
whereas if a ^ 1 this is the case of ideal current bias. 

The standard method of noise spectroscopy is to measure voltage V at the shun1j^>^. On 
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FIG. 1: Electric circuit, a) Normal tunnel junction voltage-biased through the shunt resistance R. 
b) Parallel to the shunt there is a Josephson junction, which probes the phase difference at the 
shunt, c) The IV curve of the normal junction in the limit of ideal voltage bias. 

the other hand, a Josephson junction added parallel to the shunt (Fig. ^p) may probe the 
phase difference fluctuations at the shunlj^. We assume that the both resistances, Rt and 
R, clearly exceed the quantum resistance Rk = h/e^. According to the orthodox theory^ 
at T = 0, in this limit the average rate of tunneling through the junction at Vy > e/2C is 
(Vt — e/2C)/ei?T- The IV curve for an ideal voltage bias (Vr ~ Vq) is shown in Fig. Ht. 
The junction is Coulomb blockaded as far as the voltage Vt does not exceed the Coulomb 
voltage offset e/2C. 

III. FULL STATISTICS FOR IDEAL VOLTAGE BIAS 



There is a straightforward procedure to find the full statistics in the limit of ideal voltage 
bias exactly^. Because of the constant voltage at the junction the probability of tunneling 
is also constant and the full counting statistics is given by the Poissonian distribution. The 
generating function, which yields all moments and cumulants, is 



F,(A,t) = ^e^"P„(t) = e 



(VTt/eRT)(e^-l) ^ g(/t/e)(e^-l) 



(1) 



where Pn{t) = e~^n^ /n\ is the Poissonian distribution, n is the number of tunneling events 
during the time interval t. The average number of events is determined by the current: 
ft = It/e. 

Though in the limit of ideal voltage bias the voltage at the junction V^ ~ Vq practically 
does not vary, the small voltage l^ ^ Vq at the shunt strongly fluctuates. In order to find the 



statistics of voltage V let us consider a sequence of A^ tunneling events at random moments 
of time tj {j = 1,2,...N), which may occur during the long time interval T. Time T is 
connected with A^ via the relation A^ = IT/e. Any tunneling produces the voltage pulse at 
the shunt, so the voltage varies in time as 

^ p 
V{t)=Y.re-^'-''^/^Q{t-t,), (2) 

where r = RC and G(t) is the step function. Then one can find the generating function for 
the voltage distribution averaging over random moments of time tj-. 

= ("l + Ml) V = e»l'->'/= , (3) 

where the dimensionless voltage v = VC/e was introduced, 

$,(//) = 1^" dt (e'''-''' - l) = r[Ei(-z/) - 7 - ln(-z/)] (4) 

is the contribution of a single tunneling event, 7 is Euler's constant, and Ei{z) = J^{e^ — 
l)dz/z is the exponential integral. 

Let us consider the random phase variation, which follows from Eq. Q and the Josephson 
relation dip/dt = eV/h: 

N 

^(t) = Er[l-e-(*-*^)/^]e(t-t,), (5) 

i=i 

where r = 2ttR/ Rk- Our goal is to find the full statistics for the phase difference A(/9(t) = 
^{t + to) — V'(to) = Hj 5^j during a fixed time interval t, where the contribution from the 
tunneling at the moment tj is 

% = r [1 - e-(*+*°-*^)/"] e(t + to - tj) - r [1 - g-^*""*^)/"] e(to - tj) . (6) 

Similar to the voltage statistics one should find the generating function for the phase differ- 
ence, which does not depend on to after averaging: 



^ f^dt,,.Jr (. , %{i) 



3 



N 



i^.(e,t) = (e«^^(*n = n/ ^e«'^>= 1 + :^ =,^.mi/e^ (7) 



f.Vo T 



where = yj/r is the rescaled phase and 



%{U) = -t-r{e^[E,{0-E,{^e 



-t/r 



+ Ei -ni 



-t/r 



+7 + ln[-e(l-e-*/")]} 



(8) 



We have introduced the rescaled phase since after this, statistics of phase at long times is 
identical to counting statistics (see below). Note that at ^ = ir the generating function 
F^(^,t) yields the phase correlator (^e*'^^*-'"*'^*-"-'), which determines the current through the 
normal tunnel junction in the P{E) theory, whereas at .^ = 2ir this yields the phase correlator 
for the Cooper pairs in the Josephson junction^'^. At long times t ^ t the phase statistics 
becomes identical to the counting statistics (apart from the scaling factor r): 



F^{^,t)=e 



(/i/e){e«-l) 



(9) 



This is equal to Fc{X,t), Eq. (0), with A = ^. On the other hand at short times t <^ r 
the voltage does not vary essentially during the time t, and the phase statistics should be 
identical to the voltage statistics. Indeed, at t ^ r Eq. ((7j) yields the voltage generating 
function -F^,(z/), Eq. 0, with z/ = ^t/r. 

IV. FULL VOLTAGE STATISTICS FROM MASTER EQUATION 



In order to find the full voltage statistics for an arbitrary bias we use the master equation 
for the voltage probabihty density P(V, t): 

dPiy, t) 1 d 



eRi 



^r-i.}P(V.t)-{Vr^^)p{v-lt 



dt 



RCdV 



[VP{V,t)] 







12 in the limit T 



(10) 
when 



This equation directly follows from the master equation in Ref. 

tunneling can occur only in one direction. Our analysis addresses the case of high currents 

when the voltage Vt at the junction always exceeds the Coulomb gap e/2C. We remind that 

we consider the high-impedance circuit when the resistances R and Rt exceed the quantum 

resistance. 

The equation for the generating function {y = VC/e), 



FM 



e''''P{V)dV 



(11) 



is obtained by integration of the master equation over the whole interval of relevant voltages: 



dF,{iy) ^ V dF,{v) ^ e^ - 1 



dt 



RC du 



eRi 



C du 



K)- 



2C 



FJv) 







(12) 



In terms of the dimensionless time t = t/ RC and the dimensionless voltage vq = VC/e, the 
equation becomes: 



9W^|,^„(^„_j),aF.M 



dt 



du 



a{e''-l){vo--]F,{u) 







(13) 



In the stationary case this is an ordinary differential equation with respect to u, which has 
the following solution: 



Fju) 



exp 



1 



a e 
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2 J Jo z + a(e^ — 1 



-dz 



(14) 



The solution satisfies the natural boundary condition that F^iy) = 1 at z/ = 0. According 
to Eq. (|lip this provides a proper normalization of the probability density. 

The limit a = R/Rk -^ corresponds to the ideal-bias case when the generating 
functions given by Eqs. and (fT^ are identical keeping in mind that the current is 
I = (K) — e/'^C)/{RT + R). In the general case of arbitrary a, the first three cumulants for 
the voltage at the shunt are: 

e dlnFy(u) 






_ e^ dHnF^iu) 
~C^ d^^ 



du 



IR 



IR 



2{l + a)C 
IR{2 - a) e^ 



(15) 



C3 dz/3 6(1 + a)2 C^ 

The first cumulant is simply the average voltage at the shunt. The third cumulant (skewness) 
of the voltage distribution changes sign at a = 1/2. The width of the distribution is 
determined by the smallest from the two voltages: the shunt voltage IR for the ideal voltage 
bias (a — > 0), or the junction voltage IRt for the ideal current bias (a -^ oo). 



V. FULL PHASE STATISTICS 



In order to find the full phase statistics one should consider the master equation for 
voltage and phase. We introduce the probability density PiV, ip, t) for the shunt voltage V 



and the corresponding phase (f{t) = (e/h) /* V{t')dt'. The master equation is: 

dP{V,y^,t) eVdP{V,ip,t) 1 d 



dt 



+ 



n 



dip 



1 



eRi 



Vr - i 1 P(V,^,t) -{vT + 4?^]P[V--,v,t 



2CJ 



RCdV 
2Cj V C 



[VPiV,ip,t)] 
= . 



(16) 



The full statistics of voltage and phase is determined by the generating function for the 
united phase+voltage probability distribution: 



FvA^,J^,t) 



^£,(l>+uv 



P{V,'p,t)d<pdV 



(17) 



where = ip/r is the rescaled phase and v = VC/e is the rescaled voltage. The master 
equation yields the following equation for the generating function 

^ + [u-^ + a{e - 1)] — a{e - 1) ^vq - -j F^^{^, u,t) = . (18) 

Here we use dimensionless time t = t/RC. There is a well known analogy of the phase with 
a coordinate of a diffusing particle and of the voltage with a particle velocity. Correspond- 
ingly, the phase distribution can never be stationary and permanently expands. Thus the 
generating function is always time-dependent. Performing the Laplace transformation, 



^^>^(^,^^,s) 



e-''F,^{^,iy,t)dt 



(19) 



the equation for the generating function becomes 

sF^,{^, u,s) + [u-^ + a{e^ - 1)] '^^^^^^^ ""' '^ 



di 



aie" -1) (vq- -) F^^{i, u, s) 



F,^{i,v,t)\t=^ , (20) 



where Ft,<^(^, z/, t)|t=o is the initial value of the generating function. The solution of this 
nonuniform differential equation is 



Fv,p{^,i^,s) = Fo(^, z/, s) / 



Ui+s [-^ + x + a(e^ - l)]Fo{^, x, s 
where ^i is a zero of the denominator in the integrand function. 



-dx , 



e+6 + «(e«^-i) = o. 



(21) 



(22) 



and 



Foi^,t^,s) = exp 



-s + (vo - Ij Q;(e^ - 1) 
'ii+s ~^ + x + Q;(e^ - 1) 



dx 



(23) 



is the solution of the uniform equation when the right-hand part of Eq. ()20j) vanishes. In 
order to cut divergence at the pole at a; ^ .^i, an infinitely small constant 6 was introduced. 
In final expressions this divergence is canceled and the limit 6^0 yields convergent results. 
One can check that the solution Eq. (j2ip satisfies the boundary condition Fjj^{C,, z^, t) = 1 at 
^ = 1/ = 0, which provides a proper normalization of the probability density. The Laplace 
transform of this condition is F^^{C,, z/, s) = 1/s at ^ = z/ = 0. In order to obtain the full 
phase statistics from the solution Eq. ()21|) one should choose the initial condition that 

F,^i^,u,t)\t=o = F,{iy) , (24) 

where Fy{v) is the generating function for the stationary shunt voltage distribution given 
by Eq. ()14|) . Independence of this function on ^ means that we fixed the phase y? = at 
the initial moment of time t = 0. Finally the full phase statistics after averaging over the 
voltage is given by: 

FAi, s) = F.,(f , 0, .) = F4i. 0, .) £^ ^-^——m___,, . (25) 

Taking derivatives with respect to the variable ^ conjugate to the phase difference one may 
obtain any moment or cumulant for the phase- difference probability. But even for the second 
or the third cumulant taking derivatives is a tiresome procedure in general. So we restrict 
the further analysis with some limiting cases. 

Long-time limit 

This case was also analyzed quantum- mechanically^. We shall see that our general clas- 
sical analysis completely agrees with it. 

In the case of the long-time asymptotics the most important contribution to the integral 
in pK|l comes from the close vicinity of the pole determined by Eq. ((221) • The solution of 
the uniform equation, Eq. (J2S1), can be also reduced to the contribution of the pole, and 

(so-s)/{l+a) 

I " _S' —I— _S'»-i / /) \ 

Fo(^, z/, s) ^exp 



S + So , 

; -ax 

?i+<5 -4 + a; + a(e^ - 1) 

where 

1 



,^-^ 



' - , (26) 



so={vo-^)a{e^'-l) . (27) 



For long-time asyinptotics the initial condition for the generating function is not essential, 
and one can assume for simplicity that -Fl,^(^, x, t)|t=o = F^ix) = 1 in Eq. (j^ . This means 
that at t = the voltage and phase difference at the shunt are zero. Finally the Laplace 
transform of the generating function for phase probability is approximated at long times 
with 

F^{^,s)^^—. (28) 

S — So 

After the inverse Laplace transformation the generating function in the time presentation is 

i"^(e,i)~e^o*. (29) 

This fully agrees with the result obtained by Kindermann et alA from the quantum- 
mechanical analysis. This is a manifestation of a simple law, which Kindermann et al. 
have formulated for the full statistics of two devices connected in series with the ideal 
voltage source. If the full statistics (either counting statistics or phase statistics) for any 
of two devices connected directly to the voltage source without another device are -Fi(^i) 
or -^2(^2) respectively, then for the two devices connected together the full statistics is 
F(^) = -Fi(^i) = -^2(G) at the condition that ^ = ,^1 + ,^2- In our case, the two devices are 
a normal tunnel junction and a macroscopic resistor (shunt). Then Fi{^i) = e^^*/^-"-*^ ^~^^ 
follows from the phase statistics of the junction at the ideal voltage bias, Eq. (jH)), whereas 
-^2('C2) = e*^'^*-*^^/^ corresponds to a macroscopic resistor (shunt) under fixed voltage bias IR. 
Equation (|^ is identical to the equation -Fi(^i) = F2{^ — C,i). 

The generating function Eq. ()29j) yields the following first 4 cumulants: 

d In F^ dso ~ dso d^i ~ f 1\ ~ a It 



(((p)) = r = r 1 = r 1 = r { Vn ) t = r — , (30) 

^^^" d^ d^ d^i d^ V 2/ 1 + a e ' ^ ^ 



«,.„.. .^...(,„„_^),_^...£^_J_, (3, 



^d^lnF^ _ 3/ l\~a{l-2a)_ 3/t 1 - 2a 



"'''^ - ^'^^ - ^' l"« - 2 J '1TT# - ^7(TTW • <''' 



^^^^^""^ ^e " v° 2J* {1 + ay -' e (i + «)6 • ^^^^ 
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Kindermann et al^have shown that in the hmit of ideal current bias (a ^ 1) the statistics 
of the phase at a contact of arbitrary transparency corresponds to the Pascal distributioni^. 
If the contact transparency is very low (the case of tunnel junction) the Pascal distribution 
is reduced to the chi-square distribution. One can check that this fully agrees with our 
analysis, and the found full phase statistics corresponds to the chi-square distribution for 
phase probability. 

Short-time limit 

The full phase statistics for short time intervals directly follows from a plausible as- 
sumption that the voltage does not vary during the observation time. This means that the 
generating function F^p{C,) for the phase is equal to the generating function Fy{v) for the 
voltage, Eq. ()14|) . with v = ^t/r as directly proved for an ideal voltage bias in Sec. IIIIl 

Our analysis shows that the crossover between the long-time and the short-time behavior 
is governed by the relaxation time f = CRRt/{Rt + R)- This time is different from the 
relaxation time r = RC, which characterizes the electron transport in the circuit between 
tunneling events. The two relaxation times coincide only in the ideal-voltage-bias limit. 

VI. COUNTING STATISTICS 

In order to find the full counting statistics we introduce the density Pn{V, t) of probability 
that during the time interval t n electrons tunneled through the junction and that in the 
end of the interval the voltage at the shunt is V. The master equation for this probability is 

dPn{V,t) 1 d 



+ 



dt RC dV 

1 



[VPn{V,t)] 



Vt - ^) PniV, t) - {V, + ^) ^n-1 (^ - ^^ ^)] = . (34) 



cRt 
Let us introduce the generating function: 

F,,(A, u,t) = Y. e^" / e^^Pn{V, t)dV . (35) 

The Laplace transform of the equation for the generating function is (the dimensionless time 
variable t = t/ RC was used): 

sF,„(A, z/, s) + '^^^^^y^'^ \y + a(e^+- - 1)1 - v^a^e^"^^ - 1)F,„(A, v, s) 



dv 



= F,,(A,i^,t)|i=o- (36) 
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The left-hand side of this equation is identical to the left-hand side of Eq. (PUJ) after trans- 
formation ^ = A + z/. The full counting statistics corresponds to the limit u —>■ 0. Eventually 
the generating function can be easily obtained from Eq. (j^ for the full phase statistics: 

J^i+s [-A + x + a(e^ - l)\Fo{X, x, s) 

Long times 

The only difference between the phase statistics and the full counting statistics appears 
in the initial boundary condition. Since for the long-time asymptotic the initial boundary 
condition is not essential, in this limit the cumulants for the full counting statistics can be 
obtained from those for the full phase statistics [Eqs. (J3mi3ip ] by simple rescaling. Namely, 
the relation between two types of cumulants is {{n'')) = {{(j)^)) = {{^p'^)) / r'^ . This relation 
was derived by Kindermann et al.— , but it was approximate in their quantum-mechanical 
analysis. In our purely classical approach this relation is exact. 

Short times 

In analogy with the short-time limit (t ^ f) of the phase statistics we use the fact that 
the voltage does not vary essentially during the time t. If one considers a subensemble of 
realizations, which correspond to some fixed voltage Vt at the junction during the time 
interval t, the distribution of numbers of tunneling events is Poissonian and the generating 
function is given by Eq. (^. But one should take into account that the voltage Vt at the 
junction fluctuates. Thus one should average over Vt = Vq — V, and the generating function 
is given by 

Fe(A, t) = I e(^^*/^^^)(^'-i)p(VT)ciVT = / e[(^°-^)*/^^^(^'-^)p(r)rfV . (38) 

One can see that the generating function for the counting statistics is directly connected 
with the generating function for the voltage, Eq. (fTTjl : 

F,(A, t) = e(^o)*/^^^(^'-i)F,(z>) , (39) 

where 

^ = -;^(e^-l)- (40) 
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This allows to obtain expressions for any cumulant of the full counting statistics using the 
expressions for voltage cumulants, Eq. ()15|) . The first three of them are: 



{{n)) 



d\nFc 

dX 
Vot 



Vot d In Fy dv 



eRi 



{{V)) 



cRt dv dX 

C t It 



e RtC 



(41) 
(42) 



2\\ dHuFc 



Vot d\nF„d^v rf^lnF,, ( du 



~\ 2 



eRi 



+ 



dv dX^ 

+ {{V')) 



dv"^ \dXj 



eRr 



{{V)) 



C t 
'cRtC 



RtCj 



It 

e 



a 



f 2('l 



a] 



(43) 



{{n')) 



dHnF^ 
~dP~ 






eR 



dlnF^ d^v 
du dX^ 



d'^ In F^d^udi) dHnF^ /duV 
^ di>^ dX'^dX^ dV^ \'dx] 



+3((V^))? 



t 



It 

e 



i + i 



RtCj 

3«2 



2 ^3 



t 



It 

e 

3 



e- \RtC 



a] 



t\^ a\2 

f) 6(1 + aY 



(44) 



f 2(l + a)2 

These expressions demonstrate that at short times t <^ f = CRRt/{Rt + R) the counting 
statistics becomes Poissonian even in the limit of current bias a -^ oo. Thus for short times 
the circuit (environment) effects are not so essential as for long times. It is worthwhile to 
note that "short" times t <^ f in reality are not necessarily short compared with the average 
time e/I between tunneling events since we consider the case of high currents /, and very 
many electrons may tunnel during the "short" time t <^ f. 



VII. NUMERICAL SIMULATION 

In our computational model, the time dependence of charge Q(t) = CVrit) on the June 
tion with capacitance C is obtained by integrating the equation 

dQ Vo-Q/C (dQ\ 



dt 



R 



dt 



(45) 



where the first term on the right represents charge relaxation through the shunt resistor 
i?, and the latter term represents tunneling current in the tunnel junction. According to 
Sec. iniat Vr > e/2C the average tunnehng rate (dQi/dtx) is {Q/e — l/2)/RC. Employing 
standard Monte-Carlo methods, this average rate is used to generate tunneling events, which 
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are considered to take place instantaneously on the time scales of other processes. A ready- 
made, Mersenne-Twisters-type random number generator-^ was employed in our Fortran 
code. 

The simulation was performed using parameter values close to standard mesoscopic tunnel 
junctions, namely Rt = 10 kfi and C = 1 fF which corresponds to a Coulomb voltage of 
e/2C ~ 0.1 mV. We used Vq = 5 mV for the bias voltage, while the length of the time 
record was typically set to 1 ns. With these values our time trace displayed altogether a 
few hundred tunneling events. The relaxation rate f = CRRt/{Rt + R) normally does not 
exceed 10 ps, and our simulation is basically in the long time limit discussed in Sees. El 
and IVll The time step in the integration was At = 10^^^ — 10^^'^ s. The simulation was 
initialized for 50000 iterations before starting the calculation of the actual time traces. For 
making distributions, the calculation was repeated for 5 — 10 ■ 10^ times using the previous 
simulation as an initialization for the next one. 

Fig. 2 characterizes the skewness of the distribution of the voltage V over the shunt 
resistor obtained in our simulations. The inset displays the logarithm of a voltage distri- 
bution lnP{V), calculated at a = 0.5. The asymmetry is clearly visible in the wings of 
the lnP(y)-curve. The solid curve depicts the theoretical curve for {{V"^)) obtained from 
Eq. ()15|). The master equation approach is seen to agree with the Monte-Carlo simulation 
within the scatter of the data points. 

Phase ip{t) was calculated from the simulated voltage trace by numerically evaluating 
J dtV(t) using a 3-point Simpson rule. Fig. 3 displays the result for the ratio {{f^))/{{f'^))', 
we have chosen this ratio as both the nominator and the denominator are proportional to /, 
eliminating the bias dependence off from it. The behavior at small values of a is illustrated 
in the inset in more detail. The simulated data are seen to display a change of sign in a 
similar manner as in the theory at a = 0.5. Notice that the functional dependence here 
coincides with the current fluctuation results in the low-frequency limili^. 

Fig. 4 illustrates the ratio of the skewness to the variance for the counting statistics as 
a function of the parameter a in the long-time limit, i.e. when the length of the time trace 
exceeds the relaxation time f = CRRt/{Rt + R)- In general, we find a good agreement 
between our Monte-Carlo simulation and the ratio calculated from Eqs. pip and (jH^ with 
help of the relation {{n'')) = {{ip'')) / r'^ , except for a small offset at a > 2. The inset 
magnifies the results at a < 1.5: the data on {{n^)) / {{n^)) is seen to approach the Poisson 
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FIG. 2: Skewness of the voltage distribution over the shunt resistor as a function of a = R/Rt- 
The sohd curve is calculated using Eq. (|15() . The inset displays the probability distribution PiV) 
for voltage at a = 0.5 in arbitrary units. 
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FIG. 3: Illustration of phase statistics over the shunt resistor in terms of the ratio {{(p^))/{{ip'^)). 
The solid curve is calculated using Eqs. (|31|) and H32() . The inset displays a magnification of the 
initial part of the data. 



result {{n^)) / {{n^)) = 1 as expected when a ^ 0. 



VIII. CONCLUSIONS 



Our classical approach based on the master equation provides the full statistics of electron 
transport through a normal junction in a high-impedance circuit. In addition to the full 
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FIG. 4: Ratio of the second and the third cuniulants of the counting statistics distribution. The 
sohd curve is calculated using Eqs. (|31|) and (|,32j) and the relation {{n )) = {{(p ))/r . The inset 
displays a magnification of the initial part of the data. 

counting statistics, the statistics of voltage and phase at the shunt resistor was found. The 
analysis is valid for any time interval of observation. The results are in full agreement 
with the results of the quantum-mechanical analysis performed in the long-time limit^. In 
particular, the crossover of full counting statistics from the Poissonian in the voltage-biased 
limit to the chi-square distribution in the current-biased case was obtained. The identity 
of counting statistics and phase statistics, which was found as an approximate result of 
the quantum-mechanical analysis^, was proven to be exact within classical approach. Our 
analysis shows that strong environment (circuit) effects on counting statistics become much 
weaker at short time scales (high frequencies). We have performed numerical simulations 
using Monte-Carlo methods, which fully agree with our analytical results. 
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